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ABSTRACT 


A  divergent  geostrophic  prediction  model  proposed  by  Duthie  is 
investigated  to  determine  its  usefulness  in  producing  realistic  fore- 
casts and  to  determine  the  feasibility  of  further  research  using  the 
gradient  wind. 

24-hour  forecasts  were  made  using  the  prediction  equation, 

(A)  (3)  (C)  (O) 

at  300  mb,  500  mb,  and  850  mb  for  summer  and  winter  situations,  with  and 
without  term  B  and  with  a  tuning  constant  on  term  D. 

Tables  presenting  the  error  fields  for  these  tests  are  included, 
along  with  the  plotted  forecast  and  verifying  maps.   These  tables  and 
illustrations  will  show  that  the  best  results  at  500  mb  and  300  mb  are 
obtained  with  the  B  term  out  and  a  tuning  constant  of  minus  four  and 
minus  one  respectively  on  the  D  term.  At  850  mb  different  tuning  con- 
stants produced  no  significant  changes  in  the  prognostic  charts. 

The  results  of  this  investigation  show  that  usable  prognostic  maps 
are  produced,  that  the  model  tested  has  operational  potentiality,  and 
that  further  testing  is  warranted. 
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1.    Introduction. 

With  the  advent  of  the  high  speed  digital  computer,  forecasting 
shifted  from  a  highly  subjective  art  accomplished  by  individuals  using 
hand-drawn  analyses  and  various  prognostic  methods,  many  of  which  had  no 
theoretical  foundation,  to  a  more  sophisticated  method  of  solving  the 
complicated  differential  equations  governing  the  motions  of  the  atmos- 
phere.  Many  theoretical  models  have  since  been  tested  and  the  proficien- 
cy of  the  forecaster  seems  to  have  improved.   However  many  of  these 
models  have  experienced  difficulty  in  moving  troughs  and  systems  with 
sufficient  speed. 

The  purpose  of  this  paper  is  to  introduce  a  divergent  geostrophic 
model  proposed  by  Duthie  [2].   This  model  allows  poleward  moving  parcels 
to  increase  their  relative  vorticity  and  equatorward  moving  parcels  to 
decrease  their  relative  vorticity,  thus  allowing  for  development  and  fill- 
ing.  It  will  also  move  waves  with  a  speed  of  the  basic  zonal  current. 

The  objectives  are  to  show  that  this  model,  which  uses  the  diver- 
gence of  the  geostrophic  wind,  will: 

a.  numerically  converge  and  produce  a  recognizable  weather 
map; 

b.  move  waves  at  a  speed  consistent  with  the  verifying  analyses; 
and  c.  produce  results  sufficiently  accurate  to  warrant  further 

investigation. 

The  author  wishes  to  express  his  appreciation  to  Professor  W.  D. 
Duthie  of  the  U.  S.  Naval  Postgraduate  School,  on  whose  work  this  paper 
is  based,  for  his  guidance  and  encouragement  in  this  investigation.   In 
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addition,  appreciation  is  also  expressed  to  the  personnel  of  the  U.  S. 
Fleet  Numerical  Weather  Facility  for  their  assistance  in  this  project. 
Special  thanks  are  extended  to  Mr.  Leo  Clarke  for  assistance  in  program- 
ming and  many  helpful  suggestions. 
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2.   Background. 

The  isobaric  vorticity  equation  for  frictionless  flow  is  [3] 

ifr  -  --7V-  v+~r.  (2.i) 

It  can  be  shown  [1]  that  the  divergence  of  the  geostrophic  wind 


is 


V-V?  =  -  V?  .^  U  f .  (2  2) 

Assuming 

and  substituting  the  geostrophic  divergence  for  the  divergence  in  equa- 
tion (2.1),  the  vorticity  equation  simplifies  to 


which  expanded  becomes 

DV        P        D-b  *  (2.4) 


i^_    ~    fy   4a£ 


Equation  (2.4)  shows  that  parcels  moving  poleward  will  increase 
their  relative  vorticity  and  vice  versa.  Expanding  (2.4)  further  gives 

^lh=    ~   T^'frC*.*^  -  ^T^.i^f)],         (2.5) 


where 

-I 


3*=»f""   [v^i-vz-'vUf], 


(2.6) 
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Substituting  (2.6)  into  (2.5)  gives  the  following 

(A^i  (^  lO  ^ 

Equation  (2.7)  is  a  prognostic  equation  with  a  single  variable,  z. 
Term  A  is  the  local  change,  term  B  has  the  mathematical  form  of  the 
Helmholtz  term  but  with  variable  coefficients,  term  C  is  the  advective 
term,  and  term  D  is  the  0  term  whose  sign  and  magnitude  now  depend  on 
the  relative  vorticity. 

It  is  proposed  to  test  equation  (2.7)  with  and  without  term  B  and 
with  a  tuning  constant  on  term  D. 

A  basic  linearized  expression  can  be  derived  for  equation  (2.3). 
Assuming  only  a  basic  zonal  current  U  and  that  products  of  perturba- 
tions are  small  and  can  be  ignored,  equation  (2.3)  becomes 

H'  +  u  H'  +  v'  p  =  if'  P  •  (2-8) 

Now  assuming  a  solution  of  the  form 

J  -  d0  e  1 

equation  (2.8)  reduces  to 

&  =  "U~  . 

This  shows  that  waves  will  travel  with  the  speed  of  the  basic 
current.   This  analysis  shows  that  this  model  should  move  waves  faster 
than  models  whose  waves  move  with  the  speed  of  the  Rossby  waves  [3], 

n   -  TJ  -   -^ 
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3.   Procedure. 

Data  for  850,  300,  and  500  mb  were  obtained  from  Fleet  Numerical 

Weather  Facility,  Monterey,  California,  for  summer  and  winter  situations. 

These  data  had  been  smoothed  using  a  15-  pass,  5-point  smoother  over  a 

3969-point  grid  (grid  distance  equal  to  381  km  at  60N)  using  equation  (3.1) 

S  -  D  +  A  Vs  D,  (3.i) 

where  k  is  a  constant  smoothing  factor. 

Using  this  smoothed  field,  equation  (2.7)  becomes 

which  after  applying  finite-difference  approximations  [5],  becomes 
CV2-(VZn^V)3£i  =  -y^(S,;5?)  +  ^b>n,    (3.2) 
where     "iff      ^  t  ^5  -  I'S  '  7^  ]  . 

Equation  (3.2)  was  programmed  for  the  CDC  1604  computer  using  a 
symbolic  coded  relocatable  assembly  program  (SCRAP). 

The  equation  was  solved  for  A  S  and  then  time  stepped  in  increments 
of  one  hour  using  a  forward  time  difference  (3.3)  the  first  hour  and 
centered  time  differences  (3.4)  for  succeeding  hours. 

St-t-i  =  5-t+  AS  ,  o.3) 

5-t+i  =  5t-\  \ Z  AS  ,        (3.4) 

At  the  end  of  24  hours  a  prognostic  map  was  plotted  along  with  the 
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current  analysis.   From  these  two  fields  an  error  field  was  determined 
and  from  this  error  field  a  pillow  (3.5)  and  RMSE  (3.6)were  computed  for 
all  points  north  of  the  equator. 

Pillow    =   ^-  f  A  ~  B  )  (3.5) 

N 

RMSE    =    l*liA-l»-RUow3*  (3-6) 


A  represents  the  prognostic  field,  B  the  verifying  field,  and  N  the 
number  of  points  used. 
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4.   Results. 

Table  4.1  lists  the  pillow  and  the  RMSE  for  tests  run  with  the  B 
term  of  equation  (2.7)  in  and  out  for  summer  and  winter  situations. 
All  runs  were  for  24  hours.   This  table  shows  that  the  B  term  tends  to 
increase  the  error. 

Figure  4.1  is  the  prognostic  map  with  the  B  term  in.   Superimposed 
on  this  map  in  dashed  lines  are  the  major  troughs  and  centers  of  the 
prognostic  map  with  the  B  term  out  (fig. 4. 2).  Analysis  of  this  map 
shows  that  the  B  term  does  not  affect  the  movement  of  the  waves  but 
rather  the  intensity  of  the  systems.   The  B  term  deepens  lows  and  troughs 
and  builds  highs  and  ridges.   The  high  over  Northern  Russia  is  5349m 
with  the  B  term  in  and  5321m  with  the  B  term  out.   It  verifies  as  5250m. 
The  low  over  Greenland,  Northern  Canada,  and  the  North  Pacific  covers 
a  much  broader  area  and  is  more  intense  with  the  B  term  in  than  the  same 
low  produced  with  the  B  term  out.  The  verifying  low  is  in  better  agree- 
ment with  the  low  shown  on  the  map  with  the  B  term  out. 

Another  matter  which  must  be  considered  is  the  solution  time  on  the 
computer.   With  the  B  term  in,  it  takes  105  seconds  to  solve  the  equa- 
tion for  one  hour,  90  seconds  of  which  is  used  to  relax  the  "Helmholtz" 
term  while  it  takes  only  35  seconds  to  solve  the  equation  with  the  B 
term  out,  20  seconds  of  which  is  used  in  relaxation.   Time  on  a  computer 
is  a  premium  and  must  be  considered  in  any  model. 

Figure  4.2  is  the  prognostic  map  with  the  B  term  out.   Superimposed 
on  this  map  are  the  major  troughs  and  centers  of  the  verifying  map 
(fig. 4. 3).  Analysis  shows  that  this  model  moves  the  waves  faster  than 
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B  in  vs  B  out 
500  mb 


Winter 


Summer 


Model 


RMSE 


RMSE 


B  in 
B  out 


-ll1 

1' 


266' 
212' 


0' 
-2' 


102' 
93' 


Table  4.1  RMSE  Values  for  B  term  in  vs.  B  term  out,  500  mb 
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0    00Z  08  tiEC  G5  D  -580,       '   500    MB    HEIGHTS     ,  BIN 


PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60°  NORTH  LATITUDE 
SCALE:  1:60,000.000  V- 


FLEET  NUMERICAL  WEATHER  FACILITY 
v    ,    MONTEREY,  CALIFORNIA 


Figure   4.1    24-hour  forecast  for  500  mb  from  8  December  1965,  B  term  in 
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0    80Z  08  dEC  65  D  -500.  500    MB    HEIGHTS 


PROJECTION:  POLAR  STEREOORAPHIC-TRUE  AT  60*  'NORTH  LATITUDE 
SCALE     1:60.000.000  V 


FLEET  NUMERICAL  WEATHER  FACILITY 
,        MONTEREY,  CAUrOftNlX^M^ 


Figure   4.2   24-hour  forecast  for  500  mb  from  8  December  1965,  B  term  out 
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0    00Z  09  dEC  $5  D  -500. 

PROJECTION:  POLAR  STEREOORAPHIC-TRUE  AT  60"  fcORTH  LATITUDE 
SCALE:   1:60.000,000  V- 


Figure   4.3  500-mb  analysis,  9  December  1965 


their  verifying  speeds.   Centers  also  move  too  fast  and  are  in  general 
too  intense. 

Recalling  that  in  the  perturbation  equation  (2.8),  the  p   terms 
cancelled  each  other  and  thus  resulted  in  the  waves  moving  with  a  speed 
equal  to  the  zonal  wind,  it  was  felt  that  by  changing  the  sign  of  the  y 
term  in  our  model  (D  term)  we  would  put  some  of  the  p  effect  back  in 
and  slow  the  waves  down. 

Table  4.2  lists  the  error  fields  for  summer  and  winter  situations 
at  500  mb.  This  table  shows  that  with  a  tuning  constant  of  minus  4  on 
the  D  term,  the  model  produces  the  best  results. 

Figure  4.4  is  the  prognostic  map  with  the  B  term  out  and  a  tuning 
constant  of  minus  4  on  the  D  term.   Superimposed  in  dashed  lines  are 
the  major  troughs  and  centers  of  the  verifying  map. 

Analysis  of  this  map  shows  that  the  waves  have  slowed  and  now  move 
with  a  speed  consistent  with  the  verifying  waves.   The  centers  appear 
to  be  in  quite  good  agreement  both  in  location  and  in  intensity.   It  was 
found  that  increasing  the  tuning  constant  resulted  in  continued  slowing 
of  the  waves. 

It  must  be  remembered  that  the  verifying  maps  are  prejudiced  by 
FNWF's  prognostic  model  since  this  is  used  as  a  first  guess  to  the  verify- 
ing data.  Therefore  in  sparse  data  areas  (mostly  open  oceans)  and  wher- 
ever questionable  data  are  received,  the  prognostic  data  are  used. 

Figure  4.5  is  a  500-mb  prognostic  map  for  summer  with  the  B  term 
out  and  a  tuning  constant  of  minus  4  on  the  D  term.   Figure  4.6  is  the 
verifying  map.   Of  interest  in  comparing  these  maps  is  the  agreement  on 
the  trough  coming  out  of  the  5426m  low  east  of  Greenland  and  extending 
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B  term  out 
500  mb 


Winter 


Summer 


Tuning 

8  Dec 

11 

Dec 

3 

July 

10 

July 

Constant 

P 

RMSE 

P 

RMSE 

P 

RMSE 

P 

RMSE 

1 

1' 

212' 

9' 

211' 

-2' 

93 . 

0' 

80' 

-1 

0' 

160' 

6' 

176' 

-2' 

87' 

-1' 

76' 

-2 

2' 

137' 

4' 

164' 

-2' 

84' 

-1» 

74" 

-4 

3' 

107' 

-1' 

160' 

-3' 

81' 

0' 

72' 

-8 

7' 

115' 

-4' 

90' 

30' 

140' 

Table  4.2  Results  of  RMSE  vs.  tuning  constant  for  500  mb,  B  term  out 
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8    08?  88  DEC  $5  D    500. 

PROJECTION:  POLAR  STEREOORAPHIC-TRUE  AT  60"  'NORTH  LATITUDE 
SCALE:   1:60,000,000 V 


FL&T.  NUMERICAL  WEATHER  FACILITY 
.  .._  MONTEREY.  CALIFORNIA     ,,      , 


Figure   4.4   24 -hour  forecast  for  500  mb  from  8  December  1965,  B  term  out; 
tuning  constant  minus  4 
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0    80Z  83  JUL  66  D    500.  500    MB    HEIGHTS    \  DMtN  tV 


PROJECTION:  POLAR  STEREOGRAPH1C-TRUE  AT  60'  NORTH  LATITUDE 
SCALE    1:60.000.000 ' 


FLEET  NUMERICAL  WEATHER  FACILITY 
„  .    MONTEREY,  CALIFORNIA 


chart  No.sa-t 


Figure   4.5   24 -hour  forecast  for  500  mb  from  3  July  1966,  B  term  out; 
tuning  constant  minus  4 
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0     002  04  JUL  66  D  -500 

PROJECTION:  POLAR  STEREOSRAPHIC-TRUE  AT  60*  MORTH  LATITUDE 
SCALE:   1:60.000.000  V 


FLEET  NUMERICAL  WEATHER  FACILITY 
MONTEREY,  CALIFORNIA 


Chart  No.  6B-1 


Figure   4.6   500-mb  analysis,  4  July  1966 
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down  the  European  coast  and  out  into  the  Atlantic  west  of  Gibraltar. 
Again  a  trough  in  Western  Russia  which  extends  over  the  Black  Sea  and 
into  the  Eastern  Mediterranean  compares  favorably.   The  low  over  Japan 
agrees  in  position  but  is  20m  higher  than  the  verifying  height.   The 
trough  coming  out  of  Alaska  and  down  Western  Canada  into  the  United 
States  does  not  show  a  closed  low  over  North  America  as  does  the  pro- 
gnostic map  but  a  low  of  5659m  is  indicated  on  the  verifying  map  and 
could  possibly  have  been  drawn  as  a  closed  system. 

Table  4.3  lists  the  error  fields  for  the  850-mb  level  for  both 
summer  and  winter  situations.   This  table  shows  a  substantial  decrease 
in  the  error  values.   Since  this  model  is  divergent,  it  might  be  expect- 
ed to  perform  better  at  a  level  where  divergence  has  more  effect.   It 
should  also  be  remembered  that  the  gradients  are  not  as  strong  as  at 
500mb. 

At  850mb  there  does  not  seem  to  be  any  preference  for  a  particular 
tuning  constant.   Figures  4.7,  4.8,  and  4.9  are  24-hour  forecasts  made 
from  10  July  1966  with  tuning  constants  of  minus  1,  minus  2,  and  minus 
4,  respectively.   Investigation  of  these  maps  will  show  no  significant 
difference  in  location  of  the  troughs  and  ridges  or  in  intensity  of  the 
systems.   Table  4.3  shows  that  varying  the  tuning  constant  does  not  ap- 
preciably change  the  RMSE  value  for  any  given  day. 

Figure  4.10  shows  a  typical  prognostic  map  for  850mb  and  figure 
4.11  is  the  verifying  map.   Investigation  of  these  maps  will  show  quite 
good  agreement  in  both  the  location  and  intensity  of  the  major  systems. 

Test  results,  error  fields,  and  investigations  of  the  other  maps 
in  both  summer  and  winter  show  similar  patterns  and  substantiate  these 
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B  Term  out 
850  mb 


Winter 


Summer 


Tuning 

8  Dec 

11 

Dec 

3 

July 

10 

July 

Constant 

P 

RMSE 

P 

RMSE 

P 

RMSE 

P 
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Table  4.3  Results  of  RMSE  vs.  tuning  constant  for  850  mb,  B  term  out 
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Figure   4.7   24-hour  forecast  for  850  mb  from  10  July  1966,  B  term  out; 
tuning  constant  minus  1 
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Figure   4.8   24-hour  forecast  for  850  mb  from  10  July  1966,  B  term  out; 
tuning  constant  minus  2 
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Figure  4.9   24-hour  forecast  for  850  mb  from  10  July  1966,  B  term  out; 
tuning  constant  minus  4 

31 


■■) 


t  ' 


^ 


8    00?  08  DEC  G5  D    850  850    MB    HEIGHTS        DMtN 


PROJECTION:  POLAR  STEREOGRAI'HIC -TRUE  AT  60    NOHTH  LATITUDE 
<**ai  f    i-ftnnonnnn  y 


FLEET  NUMERICAL  WEATHER  FACILITY 
.   .    MONTEREY.  CALIFORNIA 


Chart  No.  6B-1 


Figure   4.10   24-hour  forecast  for  850  mb  from  8  December  1965,  B  term  out; 
tuning  constant  minus  4 
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Figure    4.11   850-mb  analysis,  9  December  1965 
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initial  results. 

Table  4.5  is  a  comparison  between  the  error  fields  obtained  with 
this  model  and  the  error  fields  of  FNWF's  model  for  the  same  dates  at 
500  mb.   The  table  shows  that  this  model  produces  results  comparable  to 
those  of  the  present  operational  model  and  in  most  cases  might  surpass 
it  if  this  model  were  used  as  a  first  guess  field  for  analysis. 

Since  FNWF  did  not  have  an  error  field  for  850mb,  it  was  necessary 
to  look  elsewhere  for  comparison  purposes.   In  1962,  Nicholson  [4]  in- 
vestigated a  divergent  model  proposed  by  'Arnason  [1].   On  the  basis  of 
the  good  results  obtained  from  initial  testing,  •Arnason' s  model  was 
recommended  for  further  investigation.  Although  an  exact  comparison  is 
not  possible  since  the  days  chosen  for  testing  were  not  the  same,  it  is 
felt  that  some  benefit  can  be  derived  by  comparing  the  order  of  magnitude 
of  the  error  fields  for  each  model. 

Table  4.6  shows  the  850-mb  RMSE  values  for  both  summer  and  winter 
situations.   The  verification  scores  for  'Arnason' s  model  are  those 
obtained  using  the  optimum  version  of  the  model.   The  values  entered  for 
Duthie's  model  are  those  showing  the  largest  errors. 

During  initial  testing  at  300  mb,  it  was  found  that  a  relaxation 

coefficient  of  17cm  was  to  small.   This  coefficient  required  the  computer 

to  make  in  excess  of  fifty  passes  to  arrive  at  a  solution.   Since  each 

pass  takes  approximately  a  second  and  a  half,  it  required  close  to  one 

hour  to  make  a  24-hour  forecast  including  printing.   It  was  decided  to 

try  35cm  on  a  trial  basis.   This  cut  the  running  time  in  half  and  was 

used  throughout  the  remaining  tests.   However  it  is  felt  that  this 

coefficient  could  still  be  increased  and  further  experimentation  is 

recommended. 
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Table  4.4  lists  the  RMSE  values  for  tests  run  at  300  rab  for  both 
summer  and  winter  situations.   As  was  expected,  it  shows  that  the 
errors  have  increased  substantially,  being  more  than  double  those  for 
500  mb.   This  can  partially  be  explained  by  the  large  gradients  found 
at  this  level. 

Figure  4.12  is  the  24-hour  prognostic  map  from  11  December  1965  with 
the  B  term  out.   Figure  4.13  is  the  verifying  map.   Investigation  re- 
veals that  the  model  is  moving  waves  too  rapidly  eastward.   Figure  4.14 
is  the  24-hour  prognostic  map  for  the  same  time  but  with  a  tuning  con- 
stant of  minus  one  on  the  D  term.   Comparing  figure  4.14  with  figure 
4.12  shows  that  the  tuning  constant  has  slowed  the  waves  down  somewhat 
but  they  are  still  not  in  agreement  with  their  verifying  positions. 
Analysis  of  other  maps  with  various  tuning  constants  shows  that  none 
was  able  to  compensate  for  the  strong  zonal  wind  at  300  mb.   Figure 
4.15  is  the  prognostic  map  for  a  summer  situation  with  the  B  term  out 
and  a  tuning  constant  of  minus  one  on  the  D  term.   Figure  4.16  is  the 
verifying  map.   In  summer  when  the  zonal  wind  is  not  quite  as  strong, 
the  errors  are  reduced  but  the  forecast  waves  are  still  moving  too 
rapidly. 

It  is  felt  that  by  using  the  gradient  wind  instead  of  the  geostrophic 
wind  in  the  C  term  of  the  model,  it  will  reduce  the  speed  of  movement. 
Since  the  patterns  of  the  prognostic  and  verifying  maps  are  quite 
similar  this  modification  should  reduce  the  overall  error  field. 

During  the  course  of  this  investigation  it  was  found  that  the  use  of 
RMSE  to  define  an  error  field  did  not  do  an  adequate  job.  RMSE  is  only 
a  measure  of  the  difference  in  numerical  values  of  two  parameters.   It 
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Table  4.4  Results  of  RMSE  vs.  tuning  constant  for  300  mb, 
B  term  out 
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Table  4.5  RMSE  values  for  FNWF  model  vs.  Duthie  model,  500  mb 
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Table  4.6  RMSE  values  for  'Arnason  model  vs.  Duthie 
Model,  850  mb 
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Figure  4.12 


24-hour  forecast  for  300  mb  from  11  December  1965 
B  term  out  ' 
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Figure  4.13    300-mb  analysis,  12  December  1965 
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Figure  4.14   24-hour  forecast  for  300  mb  from  11  December  1965, 
B  term  out;  tuning  constant  minus  1 
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Figure  4.15   24-hour  forecast  for  300  rab  from  3  July  1966, 
B  term  out;  tuning  constant  minus  1 
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Figure  4.16   300-mb  analysis,  4  July  1966 
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gives  no  indication  as  to  how  the  maps  resemble  each  other.   It  is  felt 
that  using  RMSE  in  conjunction  with  a  correlation  coefficient,  which 
is  a  measure  of  how  much  the  products  look  alike,  would  provide  a  more 
useful  error  field.   For  example,  a  300-mb  map  compared  to  a  500-mb  map 
would  show  a  high  correlation  coefficient  but  would  also  have  a  high 
RMSE.  A  high  correlation  coefficient  together  with  a  low  RMSE  would 
indicate  that  the  prognostic  map  not  only  looked  like  the  verifying 
map  but  also  had  small  numerical  errors. 
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5.    Conclusions. 

On  the  basis  of  the  results  of  this  investigation,  the  model 
produces  a  prognostic  map  which  is  a  very  good  representation  of  the 
verifying  data.  At  both  500  mb  and  850  mb  the  waves  of  the  'tuned' 
model  move  with  a  speed  consistent  with  the  verifying  waves  and  the 
location  and  intensity  of  the  centers  agree  remarkably  well.  At  300 
mb  the  model  tends  to  move  systems  too  far.  This  is  due  to  the  use  of 
a  geostrophic  wind  with  the  strong  zonal  components  at  this  level. 

In  view  of  the  close  agreement  at  500  mb  and  below,  it  is  felt 
that  further  investigation  is  warranted,  incorporating  the  gradient 
wind  in  lieu  of  the  geostrophic  wind,  the  geostrophic  relative  vorticity 
(equation  2.6  with  the  second  term  in  the  brackets  omitted)  and  daily 
testing  begun  over  an  extended  period  to  determine  the  feasibility  of 
this  model  for  possible  operational  use. 
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